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Understanding and control of spin degrees of freedom on the surfaces of 
topological materials are the key to future applications as well as for realizing 
novel physics such as the axion electrodynamics associated with time-reversal 
symmetry breaking on the surface. We experimentally demonstrate magnet- 
ically induced spin reorientation phenomena simultaneous with a Dirac-metal 
to gapped-insulator transition on the surfaces of manganese-doped Bi2Se3 thin 
films. The resulting electronic groundstate exhibits unique hedgehog-like spin 
textures at low energies which directly demonstrates the mechanics of time- 
reversal symmetry breaking on the surface. We further show that an insulating 
gap induced by quantum tunneling between surfaces exhibits spin texture modu- 
lation at low energies but respects time-reversal invariance. These spin phenom- 
ena and the control of their Fermi surface geometrical phase first demonstrated 
in our experiments pave the way for future realization of many predicted exotic 
magnetic phenomena of topological origin. 

Since the discovery of three dimensional topological insulators [lH3, topological order 
proximity to ferromagnetism has been considered as one of the core interests of the field 
jBHTB]. Such interest is strongly motivated by the proposed time-reversal (TR) breaking 
topological physics such as quantized anomalous chiral Hall current, spin current, axion 
electrodynamics, and inverse spin-galvanic effect [914T2]. all of which critically rely on find- 
ing a way to break TR symmetry on the surface and utilize the unique TR broken spin 
texture for applications. Since quantum coherence is essential in many of these applications, 
devices need to be engineered into thin films in order to enhance or de-enhance surface-to- 
surface coupling or the quantum tunneling of the electrons. The experimental spin behavior 
of surface states under the two extreme limits, namely the doped magnetic groundstate and 
ultra-thin film quantum tunneling groundstate, is thus of central importance to the entire 
field. However, surprisingly, it is not known what happens to the spin configuration under 
these extreme conditions relevant for device fabrications. Fundamentally, TR symmetry is 
inherently connected to the Kramers' degeneracy theorem which states that when TR sym- 
metry is preserved, the electronic states at the TR invariant momenta have to remain doubly 
spin degenerate. Therefore, the establishment of TR breaking effect fundamentally requires 
measurements of electronic groundstate with a sj^m-sensitive probe. Here we utilize spin- 
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resolved angle-resolved photoemission spectroscopy to measure the momentum space spin 
configurations in systematically magnetically doped, non-magnetically doped, and ultra-thin 
quantum coherent topological insulator films [T7] , in order to understand the nature of elec- 
tronic groundstates under two extreme limits vital for magnetic topological devices. These 
measurements allow us to make definitive conclusions regarding magnetism on topological 
surfaces, and make it possible to quantitatively isolate the TR breaking effect in generating 
the surface electronic gap from many other physical or chemical changes also leading to gap- 
like behavior [THH23] often observed on the surfaces. Spin reorientation measurements and 
the systematic methodology demonstrated here can be utilized to probe quantum magnetism 
on the surfaces of other materials as well. 

Evolution of topological surface states with magnetic doping 

In order to study the evolution of topological surface states upon magnetic doping, mag- 
netically (Mn%) and (chemically similar) non-magnetically (Zn%) doped Bi2Se3 thin films 
are prepared in high quality using the molecular beam epitaxy (MBE) growth method. 
A sample layout and a photograph image of a representative MBE grown film used for 
experiments are shown in Fig. la and b. Using standard X-ray magnetic circular dichro- 
ism [2U [23], we characterize the magnetic properties of doped Bi2Se3 films (Fig. Ic). In 
Mn-doped Bi2Se3, a hysteretic behavior in the out-of-plane magnetic response is observed 
consistently which suggests a ferromagnetically ordered groundstate. The observation of 
ferromagnetic character and its absence in Zn-Bi2Se3 motivate us to systematically compare 
and contrast the electronic density of state behavior in the vicinity of the surface Dirac node 
of these samples. Fig. Id shows the measured electronic states of Mn(Zn)-doped Bi2Se3 us- 
ing high-resolution (spin-integrated) angle-resolved photoemission spectroscopy (ARPES). 
In the undoped Bi2Se3 film (leftmost panel of Fig. Id), a map of spectral density of states 
reveals a bright and intact Dirac node (signaled by the red spot located at the Dirac crossing 
in the panel), which indicates that in undoped Bi2Se3 the Dirac node is gapless, consistent 
with the previous studies of pure Bi2Se3 [26]. In samples where Mn atoms are doped into the 
bulk (first row of Fig. Id), we observe that the corresponding bright (red) spot at the Dirac 
node gradually disappears, revealing a clear systematic spectral weight suppression (SWS) 
with increasing Mn concentration. In contrast, the spectral intensity at the Dirac node is 
observed to survive upon systematic Zn doping except for the Zn = 10% sample where some 
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suppression of intensity is observed. This suggests that the Dirac node remains largely intact 
upon Zn doping. The observed systematic behavior of spectral evolution motivates us to 
quantitatively define an energy scale, Esws; associated with the SWS observed at the Dirac 
node. The Esws is taken as the energy spacing between the upper Dirac band minimum 
and the Dirac node location along the energy axis as illustrated in Fig. le, which roughly 
corresponds to half of the surface gap magnitude. The value of the energy scale can be 
quantitatively determined by fitting the ARPES measured energy- momentum distribution 
curves (method described in supplementary section VI) . The doping dependence of the Egws 
on samples measured at T = 20 K is shown in Fig. 2c. The Esws is observed to increase 
rapidly with Mn concentration but it remains nearly zero with Zn doping. The effect of 
temperature dependence on Esws is shown in Fig. 2d. The temperature induced decrease 
of Esws is consistent with gradual weakening of magnetism. These observations collectively 
reveal direct correlation between X-ray magnetic circular dichroism measured ferromagnetic 
character and the ARPES measured SWS (or gap) on the Mn-Bi2Se3 films. 

Although the predominant trend in the doping evolution of surface states suggests a 
correlation between magnetism and SWS, we do notice a non- negligible Esws at high Zn 
concentration (10%). This is likely due to increasing chemical disorder on the surface of 
the film since disorder degrades the surface quality. Similarly, magnetization measurements 
show that ferromagnetism vanishes before reaching T = 300 K (inset of Fig. 2b), but nonzero 
SWS (Fig. 2d) is still observed perhaps, in this case, due to thermal disorder of the relaxed 
film surface at high temperatures. Therefore, the correlation between ferromagnetism and 
the ARPES gap is only clear at low temperatures in samples with reduced disorder. The 
momentum width (Ak) of the surface electronic states can be taken as a rough relative 
measure of surface disorder, sample to sample, which is found to significantly increase upon 
both magnetic and non- magnetic dopings (Fig. 2c). Strong spatial fluctuations of the surface 
electronic states in doped Bi2Se3 has been observed in a recent STM work by Beidenkopf 
et al. [18], where the authors suggest the observation of gap-like feature at the Dirac point 
without breaking TR symmetry. These ambiguities |2T] associated with the observed gap-like 
feature across many different experiments strongly call for critically important spin-resolved 
measurements which also serve as a collective method, as we show, to unambiguously identify 
the correct nature of the gap. 
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Spin configuration of the magnetic groundstate 

In order to study the evolution of spin texture upon magnetic doping, we perform spin- 
resolved measurements on Mn-Bi2Se3 topological surface states. We present two indepen- 
dent but representative spin-resolved ARPES measurements on Mn(2.5%)-Bi2Se3 film I and 
film IL Films I and II, both containing same nominal Mn concentration, are measured and 
analyzed using two different spectroscopic modes, namely, spin-resolved momentum distri- 
bution curve (spin-resolved MDC) measurement mode and spin-resolved energy distribution 
curve (spin-resolved EDC) measurement mode (see Methods section), in order to exclude 
any potential systematic error in the spin measurements. Fig. 3a-d show measurements 
on film I. Our data shows that out-of-plane spin polarization is nearly zero at large 
momentum k// far away from the Dirac point energy (0 < Eb < 0.1 eV in Fig. 3c, d). 
While approaching the Dirac point (0.1 eV< Eb < 0.3 eV), an imbalance between the spin- 
resolved intensity in +z and —z is observed (Fig. 3c). The imbalance is found to become 
systematically more pronounced in the data set where scans are taken by lowering the en- 
ergy toward the Dirac point. This systematic behavior observed in the data reveals a net 
significant out-of-plane spin polarization in the vicinity of the "gapped" Dirac point or near 
the bottom of the surface state conduction band. More importantly, the out-of-plane spin 
component P^ does not reverse its sign in traversing from —kji to Such behavior 

is in sharp contrast to the spin textures observed in pure Bi2Se3 [28] where spins point to 
opposite directions on opposite sides of the Fermi surface as expected from TR symmetry. 
Therefore, our P^ measurements on film I near the gap edge reveal the TR breaking nature 
of the Mn-Bi2Se3 sample where magnetic hysteresis was observed using X-ray magnetic cir- 
cular dichroism technique. In order to directly measure the spin of the surface state at P 
(the Kramers' momentum, kjj = 0), we perform spin measurements on Mn-Bi2Se3 film II 
(same Mn concentration as film I) working in the spin-resolved EDC mode. The measured 
out-of-plane spin polarization (P^) is shown in Fig. 3g. We focus on the P^ measurement 
at P, the Kramers' momentum k// = (red curve): the surface electrons at TR invariant 
P are clearly observed to be spin polarized in the out-of-plane direction. The opposite sign 
of P^ for the upper and lower Dirac band (red curve in Fig. 3g) shows that the Dirac 
point spin degeneracy is indeed lifted up (E(A;// = 0,t)7^E(^// = 0,|)). Such observation 
directly counters the Kramers' theorem and therefore manifestly breaks the TR symmetry 
on the surface. Next we analyze P^ measurements at finite k// (green curves in Fig. 3g) 
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to extract the detailed configuration of tlie spin texture. In going to larger fc// away from 
the r momenta, the measured is found to gradually decrease to zero. Moreover, the 
constant energy momentum space plane at the Dirac point (Eb = Ed) is observed to serve 
as a mirror plane that reflects all of the out-of-plane spin components between the upper 
and lower Dirac bands. Thus both spin-resolved MDC (film I) and spin-resolved EDC (film 
II) measurement modes result in the consistent conclusions regarding the spin configuration 
of the films. These systematic measurements, especially at the vicinity of the gap, reveal 
a hedgehog-like spin configuration for each upper (or lower) Dirac band separated by the 
magnetic gap, which breaks TR symmetry, as schematically presented in the inset of Fig. 
3f (see supplementary for additional data and analysis). 

Spin configurations of the non-magnetic groundstates 

Spin texture measurements on non-magnetically doped films Zn(1.5%)-Bi2Se3 are presented 
in Fig. 4a-d. The out-of-plane polarization measurements reveal a sharp contrast to the 
magnetically doped Mn-Bi2Se3 films, specifically, the near absence of finite P^ component 
around T within our experimental resolution (Fig. 4d). A very small P^, however, at large fc// 
is observed, which is expected due to surface state warping also observed in other topological 
insulator (TI) compounds |17] (Fig. 4d). The signal being associated with warping is further 
confirmed in our data due to their TR symmetric nature, that P^ is observed to reverse its 
sign in traversing from —k// to +k//. Moreover, our in-plane spin measurements (Fig. 4c) 
show that Zn-Bi2Se3 film retains the helical spin texture protected by the TR symmetry, 
as observed in pure Bi2Se3 and Bi2Te3 single crystals [28j. Therefore we conclude that 
nonmagnetic Zn doping does not induce observable spin reorientation on the topological 
surface. The contrasting behavior observed between Mn-Bi2Se3 and Zn-Bi2Se3 samples as 
presented in Fig. 3 and 4 provide clear evidence for TR symmetry breaking in Mn-Bi2Se3. 

A surface band gap at the Dirac point can also be generated in Bi2Se3 in its ultra-thin 
film limit. In this case, the top and bottom surfaces couple together and open up a gap as 
electrons can tunnel from one to the other (Fig. 4e-h). Such a gap in the surfaces is not 
related to magnetism. It is important to know the spin configuration associated with such 
a gap. In Fig. 4e-h, we utilize Spin-resolved ARPES to measure the spin configuration on 
the very top region (within 5 A) of a Bi2Se3 film whose thickness is three quintuple layers (3 
QL 28.6 A). At large parallel-momenta far away from f (e.g. -0.10 A ^ in Fig. 4g), we 
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observe clear spin polarization following left-handed helical configuration with the magnitude 
of the polarization around 35 — 40%. However, in going to smaller fc//, the magnitude of the 
spin polarization is observed to be reduced. At the TR invariant T momenta, spin-resolved 
measurements (Fig. 4g red curve) show no net spin polarization. This reduction of the spin 
polarization at small momenta near the gap is an intrinsic effect (see supplementary Fig. 
S19). These observations can be understood by considering the scenario where the surface- 
to-surface coupling dominates, and the two energetically degenerate surface states from top 
and bottom that possess opposite helicities of the spin texture cancel each other at T [T7] . 
This results in strong suppression of spin polarization in the vicinity of this gap, whereas 
upon probing momenta to large fc// away from f, the finite kinetic energy of the surface 
states (oc vk//) naturally leads to the spatial decoupling of two Dirac cones. These spin 
measurements on the ultra-thin Bi2Se3 film reveal the interplay between quantum tunneling 
(coupling) and the spin texture modification, which is of importance in spin-based device 
design with thin films. The observed spin texture however does not break TR symmetry, 
since the spins remain doubly degenerate at the TR invariant momenta T. This is in clear 
contrast to the spin texture observed in Mn-Bi2Se3. 

Magnetic contribution and geometrical phase tuning 

The magnetic contribution to the gap of the Mn-Bi2Se3 film can be quantitatively identified 
using the spin texture data. The simplest k-p Hamiltonian that describes topological surface 
states under TR symmetry breaking can be written as H = v{kxay — kyax)+bzO'z, where a and 
k are the spin and momentum operators respectively, bz corresponds to half of the magnetic 
gap and v is the velocity of the surface Dirac band. We specify the out-of-plane polar angle 
9 of the spin polarization vector (inset of Fig. 3g) as 6' = arctan^. The magnitude of 
the polar angle 6 reflects the competition between the out-of-plane TR breaking texture (oc 
bz) and the in-plane helical configuration component (oc vk//). Using the measured spin- 
resolved data sets {6, k), we fit the magnetic interaction strength bz within a k-p scenario (see 
supplementary section II. 2). As an example, we fit bz based on spin-resolved data sets in Fig. 
3g on Mn(2.5%)-Bi2Se3 film II, as shown in Fig. 5c and obtain a value of 21 meV. This is a 
significant fraction of the SWS energy scale observed on the same sample, Esws > 50 meV 
(see Fig. 2c for Mn = 2.5%) obtained from the spin-integrated measurements in Fig. Ic. 
Thus we identify the magnitude of the magnetic contribution (bz) to the observed spectral 
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weight suppression using spin-sensitive measurements, which suggests that the magnetic 
contribution is significant to Esws in Mn-Bi2Se3. 

As demonstrated recently [22], the geometrical phase (GP) defined on the spin texture 
of the surface state Fermi surface [3] (also known as the Berry's phase) bears a direct cor- 
respondence to the bulk topological invariant realized in the bulk electronic band structure 
via electronic band inversion |28l [29]. We experimentally show that a GP tunability can 
be realized on our magnetic films which is important to prepare the sample condition to 
the axion electrodynamics limit. On the Mn-Bi2Se3 film, spin configuration pattern can be 
understood as a competition between the out-of-plane TR breaking component and the in- 
plane helical component of spin. The in-plane spin that can be thought of winding around 
the Fermi surface in a helical pattern contributes to a nonzero GP [28] , whereas the out-of- 
plane TR breaking spin direction is constant as one loops around the Fermi surface hence 
does not contribute to the Berry's phase (GP). As a result, the GP remains almost vr if 
the chemical potential lies at energies far away from the Dirac point, whereas it starts to 
decrease and eventually reach to as one approaches the TR breaking gap by lowering the 
chemical potential as discussed in theory [30], at least within the magnetic energy scale bz 
(Fig. 5). We show that this theoretical requirement can be experimentally achieved on the 
Mn-Bi2Se3 surface via surface NO2 adsorption at some specific dosage level. Fig. 5a shows 
the Mn(2.5%)-Bi2Se3 surface states with in situ NO2 adsorption. The chemical potential 
is observed to be gradually shifted and finally placed within the magnetic gap. The asso- 
ciated phase (GP) at each experimentally achieved sample chemical potential is found to 
gradually change from vr to 0. The GP = is the experimental condition for realizing axion 
electrodynamics with a topological insulator [9l [TT] . 

With the chemical potential moved into the magnetic gap, the time-reversal breaking in- 
gap state features a singular hedgehog-hke spin texture (Fig. 5d). Such spin configuration 
simultaneous with the chemical potential placed within the magnetic gap (Fig. 5d) is the 
fundamental requirement for most of the theoretical proposals relevant to the utilization of 
magnetic topological insulators in novel devices [H [9HT2] . Furthermore, if the bulk band- 
gap of the Mn-Bi2Se3 film is tuned to zero at the critical point of the topological phase 
transition [29] , a new topologically protected Weyl semimetal phase [311 [32| with yet more 
exotic but undiscovered state of matter is also predicted to take place, which is among the 
most exciting future frontiers to be enabled by our achievement of a sample that features a 
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TR broken hedgehog-like spin texture with a GP = state [311 132] . 
Methods 

Sample growth: The Mn doped Bi2Se3 thin films were synthesized by Molecular Beam 
Epitaxy using high purity elemental (5N) Mn, Bi, and Se sources. A thin GaAs buffer layer 
was first deposited on the epi-ready GaAs 11 lA substrate after thermal desorption of the 
native oxide under As pressure. Then the substrate was transferred to another chamber 
without breaking the vacuum, where a second buffer layer of ZnSe was deposited for further 
smoothing the surface. Mn-Bi2Se3 layer (~ 60 nm) was then grown with a high Se/Bi beam 
equivalent pressure (BEP) ratio of ~ 15. The Mn doping concentration was controlled 
by adjusting the Bi/Mn BEP ratio ranging from 8 to 60. To protect the surface from 
oxidation, a thick Se capping layer was deposited on the Mn-Bi2Se3 thin film immediately 
after the growth. The Zn doped Bi2Se3 control samples were also synthesized under the 
same conditions as Mn-Bi2Se3, with the Zn doping concentration controlled by Bi/Zn BEP. 
Ultra-thin Bi2Se3 film was prepared in thickness of 3 QL, with a typical 1 QL peak to peak 
variation (2 — 4 QL). Complete details of film growth are presented in Ref [33] and also 
detailed in the supplementary. 

Electronic structure measurements: Angle-resolved photoemission spectroscopy 
(ARPES) measurements (spin-integrated) were performed with 29 eV to 64 eV photon en- 
ergy on beamlines 10.0.1 and 12.0.1 at the Advanced Light Source (ALS) in Lawrence Berke- 
ley National Laboratory (LBNL). Spin-resolved ARPES measurements were performed on 
the 13 beamline at Maxlab in Lund, Sweden, COPHEE spectrometer SIS beamline at the 
Swiss Light Source (SLS) in Switzerland, and UE112-PGM1 beamline PHOENEXS chamber 
at Bessyll in Berlin, Germany, using photon energies of 8 — 11 eV, 20 — 22 eV, and 55 — 60 
eV for the three beamlines respectively. For ARPES measurements, samples were annealed 
in situ to evaporate the amorphous Se cap layer in order to reveal the clean surface (see 
supplementary for further details). 

Spin-resolved measurements: Spin-resolved measurements were performed using 
spin-resolved ARPES beamlines mentioned above with double classical Mott detectors and 
linearly polarized photons. The spin-resolved measurements in Fig. 3 and 4 were per- 
formed in two different modes: spin-resolved MDC mode and spin-resolved EDC mode, 
spin-resolved MDC mode means each single measurement measures the spin polarization at 
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a fixed binding energy (Eb) along certain momentum cut direction in momentum space (see 
Fig. 3a-d on Mn(2.5%)-Bi2Se3 film I), spin-resolved EDC mode means each single measure- 
ment measures the spin polarization at a fix momentum (k) along the binding energy axis 
(Eb) (see Fig. 4e-g and Fig. 4 on Mn(2.5%)-Bi2Se3 film II, Zn(1.5%)-Bi2Se3 film, and 3 
QL undoped Bi2Se3 film). The details of the spin-resolved technique and principles can be 
found in supplementary. 

Magnetic property characterizations: In order to cross-check the ferromagnetism 
implied by the spin measurements on Mn-Bi2Se3 surface. X-ray magnetic circular dichroism 
(XMCD) was independently performed on Mn-Bi2Se3 film surface at the back-endstation 
of DlOll beamline at Maxlab in Lund, Sweden, with total electron yield (TEY) mode at 
temperatures ranging from 40 K to 300 K. XMCD measurements were performed on the 
L23 absorption edge of the Mn atom by the standard methods [21], which is widely used to 
study the magnetic properties of transition metals and dilute magnetic semiconductor thin 
films or monolayers [23] . 

The detailed principles and methods of all techniques used in the experiments including 
sample growth, ARPES and spin-resolved ARPES, NO2 surface adsorption, and XMCD are 
further elaborated in the supplementary information. 
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FIG. 1: Magnetic (Mn%) and non-magnetic (Zn%) doping on Bi2Se3 films, a and 

b, Schematic layout and photograph of MBE grown Bi2Se3 films used for the experiments, c, 
Magnetization measurement at T = 45 K using magnetic circular dichroism shows out-of-plane 
ferromagnetic character of the Mn-Bi2Se3 film (111) surface through the observed hysteretic re- 
sponse (detailed in supplementary section II. 5. 2). Inset shows the measurement geometry. L(R)CP 
represents left(right)-handed circularly polarized light, d, Electronic band dispersion of Mn(Zn)- 
doped Bi2Se3 MBE thin films along the M — f — M momentum space cut. The doping level (noted 
in the top-right corner of each panel) reflects nominal Mn(Zn) concentration, which is defined as 
the ratio of Mn(Zr^+Bi o'^^r the entire film crystal. Increasing Mn concentration leads to spectral 
weigh suppression at the Dirac point, signaling a "gap". The "gap" and the hysteresis are possibly 
correlated in the data, e, Energy-momentum distribution curves of Mn(Zn)-doped Bi2Se3 samples. 
The energy scale associated with the spectral weight suppression (SWS) Esws is observed as the 
energy spacing between the upper Dirac band minimum and the Dirac point location along the 
energy axis. 
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FIG. 2: Temperature and doping dependence of magnetically induced changes on Mn- 
Bi2Se3 surface, a, The Mn atoms on the surface of the film are out-of-plane magnetically ordered, 
serving as a local magnetic field which results in the spin texture (see Fig. 3) reorientation, b, Two 
independent hysteresis measurements at T = 45 K using X-ray magnetic circular dichroism reveal 
the ferromagnetic character of the Mn-Bi2Se3 film surface. The lower inset shows the remanent 
surface magnetization as a function of temperature. The out-of-plane magnetic hysteresis and 
ARPES gap were found to be correlated with each other. The upper inset shows the gap at 
the Dirac point in Mn(2.5%)-Bi2Se3 film, c, The spectral weight suppression energy scale Esws 
(defined in Fig. 1) and inverse momentum width 1/Ak of the surface states are shown as a function 
of Mn and Zn concentration measured at T = 20 K. d, Temperature dependence of spectral weight 
suppression energy scale around the Dirac point of Mn(2.5%)-Bi2Se3 film (as noted in panel c by 
the dotted square). The Esws decreases as temperature is raised signaling gradual weakening of 
magnetism. 
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FIG. 3: Spin configuration measurements of a magnetic topological insulator. Panels a- 
d present spin-resolved measurements on film I using 20 eV photons in the momentum distribution 
curve (MDC) mode. Panels e-g present measurements on film II using 9 eV photons in the energy 
distribution curve (EDC) mode, a and b, Spin-integrated data and corresponding MDCs. c, 
Spin-resolved MDC spectra for out-of-plane direction as a function of electron binding energy, d, 
Measured out-of-plane component of the spin polarization, presented in terms of respective out- 
of-plane polar angles (6) defined in the inset of g. e-f, Spin-integrated dispersion and EDCs. 
The EDCs selected for detailed spin-resolved measurements are highlighted in green and red. g, 
Measured out-of-plane spin polarization of the EDCs corresponding to panel f. Inset defines the 
definition of spin polarization vector P and the polar angle 9. The momentum value of each 
spin-resolved EDC is noted on the top. The polar angles (9) of the spin vectors obtained from 
measurements are also noted. The 90° polar angle observed at T suggests that the spin vector at T 
momenta points in the vertical direction. The spin behavior at T and its surrounding momentum 
space reveals a hedgehog- like spin configuration for each Dirac band separated by the gap, which 
breaks time-reversal symmetry (E{k = 0,t)7^E(A; = 0,4,)), as schematically shown in the inset of f. 
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FIG. 4: Spin configurations on non-magnetic samples, a-d, Spin-resolved measurements 
on 1.5% non-magnetic Zn-Bi2Se3 film. The in-plane polarization measurements (c) reveal the 
helical spin configuration, as in pure Bi2Se3 topological insulator [28], suggesting that non-magnetic 
impurities do not induce spin reorientation on the topological surfaces. Out-of-plane measurements 
(d) show that no significant out-of-plane spin polarization is induced near the T point (the time- 
reversal invariant momenta), leaving the system time-reversal invariant overall, e-g, Spin-resolved 
ARPES measurements on ultra-thin undoped 612863 film of three quintuple- layer thickness. The 
net spin polarization is found to be significantly reduced near the gap edge around the T momenta. 
This is consistent with the fact that in ultra-thin films electrons tunnel between the top and bottom 
surfaces, h and i, A schematic of the two types of spin textures observed in our data. 
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FIG. 5: Chemical potential tuned to lie inside the magnetic gap. a, Measured surface 
state dispersion upon in situ NO2 surface adsorption on the Mn-Bi2Se3 surface. The NO2 dosage 
in the unit of Lanemuir (IL = 1 X 10"*^ torr •sec) and the tunable geometrical phase (see text) 
associated with the topological surface state are noted on the top-left and top-right corners of the 
panels, respectively. The red arrows depict the time-reversal breaking out-of-plane spin texture at 
the gap edge based on the experimental data, b, Geometrical phase (GP) associated with the spin 
texture on the iso-energetic contours on the Mn-Bi2Se3 surface as a function of effective gating 
voltage induced by NO2 surface adsorption. Red squares represent the GP experimentally realized 
by NO2 surface adsorption, as shown in a. GP = (NO2 = 2.0 L) is the condition for axion 
dynamics, c, The magnetic interaction strength bz (see text for definition), which corresponds to 
half of the magnetic gap magnitude, is obtained based on spin-resolved data sets (polar angle 9, 
momentum k) for Mn(2.5%)-Bi2Se3 film II (see Fig. 3e-g). d, The time-reversal breaking spin 
texture features a singular hedgehog-like configuration when the chemical potential is tuned to lie 
within the magnetic gap, corresponding to the experimental condition presented in the last panel 
in a. e and f, Spin texture schematic based on measurements of Zn-doped Bi2Se3 film (60 QL), 
and 3 QL undoped ultra-thin film with chemical potential tuned onto the Dirac point energy or 
within the tunneling gap. 
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SI I. Materials and methods 

Spin- integrated angle-resolved photoemission spectroscopy (ARPES) measurements were 
performed with 29 eV to 64 eV photon energy on beamlines 10.01 and 12.01 at the Advanced 
Light Source (ALS) in Lawrence Berkeley National Laboratory (LBNL). Spin-resolved 
ARPES measurements were performed on the 13 beamline at Maxlab [Ml [35 ] in Lund, Swe- 
den, COPHEE spectrometer SIS beamline at the Swiss Light Source (SLS) in Switzerland 
[361 EI], and UE112-PGM1 beamline PHOENEXS chamber at Bessyll in Berlin, Germany, 
using the classical Mott detectors and photon energies of 8 — 11 eV, 20 — 22 eV, and 55 — 60 
eV for the three beamlines respectively. For ARPES measurements, samples were decapped 
in situ to prepare the clean surface (see below), and measured at temperature below 50 K 
at chamber pressure less than 2x10^^° torr at all beamlines. 

The Mn doped Bi2Se3 thin films were synthesized by Molecular Beam Epitaxy (MBE) 
using high purity elemental (5N) Mn, Bi, and Se sources. A thin GaAs buffer layer was 
first deposited on the epi-ready GaAs lllA substrate after thermal desorption of the native 
oxide under As pressure. Then the substrate was transferred to another chamber without 
breaking the vacuum, where a second buffer layer of ZnSe was deposited to further smooth 
the surface. Mn-Bi2Se3 layer (~ 60 nm) was then grown with a high Se/Bi beam equivalent 
pressure (BEP) ratio of ~ 15. The Mn doping concentration was controlled by adjusting 
the Bi/Mn BEP ratio ranging from 8 to 60. To protect the surface from oxidation, a thick 
Se capping layer was deposited on the Mn-Bi2Se3 thin film immediately after the growth 
(see Fig. [T^). The Zn doped Bi2Se3 control samples were also synthesized under the 
same conditions as Mn-Bi2Se3, with the Zn doping concentration controlled by Bi/Zn BEP. 
Ultra-thin Bi2Se3 film was prepared in thickness of 3 QL, with a typical 1 QL peak to peak 
variation (2 ~ 4 QL). 

In order to reveal the clean Bi2Se3 surface needed for photoemission measurements, the 
MBE-grown thin films are heated up inside the ARPES chamber to ~ 250°C under vacuum 
better than 1 x 10~^ torr to blow off the Se capping layer on top of the Bi2Se3 film. Fig. 
demonstrates the decapping process by ARPES core level spectroscopies measurement. Only 
selenium core level is observed before the decapping process (blue curve), whereas both sele- 
nium and bismuth peaks are shown after the decapping (red curve). ARPES measurements 
are then performed on the clean thin film surface. Fig. ^-d show a typical spin-integrated 
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Supplementary Figure 1: Demonstration of sample surface preparation (decapping) 
procedure in ARPES measurements, a, Sample layout of the MBE grown Mn(Zn)-doped 
Bi2Se3 films, b, Core level spectroscopies on MBE thin film before and after the decapping proce- 
dure, c-d, typical spin-integrated ARPES dispersion mapping and corresponding energy dispersion 
curves (EDCs) of undoped and ~ 1% Mn-doped (nominal concentration) Bi2Se3 thin film along 
high symmetry M — F — M direction. 

ARPES measured dispersion mappings after decapping of both undoped Bi2Se3 and ~ 1% 
(nominal concentration, same as the maintext) Mn-doped Bi2Se3 thin film. 

SR-ARPES measurements are performed with double Mott detector setup, which enables 
to measure spin-resolved photoemission spectra along momentum fc// at fixed binding energy 
Eb (SR-MDC) or along binding energy Eb at fixed momentum kjj (SR-EDC) for all three 
spatial directions (Kfi) -^yti) ^"^^ -^-^ti)- The detailed and comprehensive review of SR- 
ARPES technique is presented in Ref. [3H]- The spin-resolved data on Mn-Bi2Se3 shown in 
the main paper is done in remanent magnetization mode: the film is pre-magnetized in the 
UHV preparation chamber in the ARPES system at the base temperature by approaching 
the film with a permanent magnet with magnetic field roughly at 0.4T, in order to induce 
the remanent magnetization. The film is then transferred from the preparation chamber to 
the measurement chamber to perform spin-resolved measurements, maintaining the pressure 
and temperature. 

The surface sensitivity of ARPES measurements can be obtained from the photoemission 
electron escape depth. The electron escape depth with incident photon energy from (50 — 60 
eV) are well below 5 A. This is important to explain the SR-ARPES data on the ultra-thin 
3 QL Bi2Se3 in the maintext, because a < 5 A surface sensitivity means SR-ARPES probes 
the very top quintuple layer of the 3 QL film, whose total thickness is around 3 nm. 

Magnetization of the surface in the presence of Mn atoms are checked with X-ray magnetic 
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circular dichroism (XMCD) using total electron yield (TEY) mode that greatly enhances 
the surface sensitivity (within 50 A) of the magnetization measurements. XMCD measures 
the dichroism signal of the L23 absorption edge spectra, which is widely used to study the 
magnetic properties of transition metals and dilute magnetic semiconductor thin films or 
monolayers. The X-ray magnetic circular dichroism (XMCD) measurements are performed 
at the back-endstation of DlOll beamline at Maxlab in Lund, Sweden with total electron 
yield (TEY) mode on Mn-Bi2Se3 film. The temperature of the XMCD measurements ranges 
from 40 K to 300 K. The magnetic field is applied in the out-of-plane direction normal to the 
sample surface, with a maximum H field of 300 oe. The left (right) hand circularly polarized 
incident light is also normal to the film surface, with an accuracy of the circular polarization 
better than 90%. The sample is prepared the same way as in ARPES measurements in 
situ by decapping inside the XMCD UHV chamber. The XMCD technique and the XMCD 
measurements on Mn-Bi2Se3 film surface are discussed in detail in the section below SI. II. 
3. 

Adsorption of NO2 molecules on sample surface was achieved via controlled exposures to 
NO2 gas (Matheson, 99.5%). The adsorption effects were studied under static flow mode by 
exposing the clean sample surface to the gas for a certain time at the pressure of 1 x 10^^ 
torr, then taking data after the chamber was pumped down to the base pressure. Spectra 
of the NO2 adsorbed surfaces were taken within minutes of opening the photon shutter to 
minimize potential photon induced charge transfer and desorption effects. 

The first-principles calculations are based on the generalized gradient approximation 
(GGA) [39] using the full-potential projected augmented wave method as implemented 
in the VASP package [H]. A 2 x 2 x 6 quintuple-layer slab model with a vacuum thickness 
lager than 10 A is used in this work. To simulate the dilute Mn doping, one Bi atom is 
substituted by one Mn atom in both the top and bottom layers with the internal atomic 
coordinates optimized. The electronic structure calculations were performed over 5x5x1 
Monkhorst-Pack k-mesh with the spin-orbit coupling included self-consistently. 
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SI II. Mn-Bi2Se3 



1. Comprenhensive spin-resolved measurements 

a Mn-Bi,Se, film I 
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Supplementary Figure 2: In-plane spin measurements of Mn-Bi2Se3 film I. a, Spin- 
integrated ARPES dispersion mapping at incident photon energy 20 eV (used for spin-resolved 
measurements), b-f, In-plane spin polarization along x and y directions, g, The x direction is 
defined as along the spin- resolved MDC measurement direction. 



Here we show comprehensive spin-resolved data on Mn-Bi2Se3 film I and II in the 
maintext Fig. 3. Mn-Bi2Se3 film I is measured using the COPHEE spectrometer 
SIS beamhne at the Swiss Light Source (SLS) in Switzerland (Fig. [2|i) [361 EZ]- The 
measurements are performed in spin-resolved momentum distribution curve (MDC) 
mode (along momentum, with fixed energies, as shown in Fig. |2^). Since the out-of- 
plane component is already presented in the maintext (Fig. 3a-d), here we show the 
in-plane components of the spin-resolved measurements. 

The X direction is defined as along the spin-resolved MDC measurement direction. 
Therefore, if the SR-MDC measurements perfectly cross the Brillouin zone (BZ) center 
r point, then polarization along x direction should be strictly zero, whereas Py 
should be opposite at the opposite side of the Fermi surface representing the helical 
spin texture [28]. Fig. |2b-f show the measured in-plane spin polarization along both x 
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and y directions. Indeed, Py reverses its sign when going from —k to +k. Meanwhile, 
a finite P^^ is also observed which becomes increasingly significant when going to large 
binding energies close to the Dirac point. Such nonzero P^^ is owning to a small tilt of 
the sample, which results in a "miss-cut" away from the BZ center P point (see Fig. [2^). 
The observed P^^ (opposite sign at opposite sides of the Fermi surface) and P^, (same 
sign at opposite sides of the Fermi surface) proves that the in-plane components of the 
spin texture still preserve the helical spin configuration. 
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Supplementary Figure 3: In-plane spin measurements of Mn-Bi2Se3 film II. a-b, Spin- 
integrated ARPES measured dispersion and EDCs with incident photon energy 9 eV (used for 
spin-resolved measurements). The EDCs selected for spin-resolved measurements are highlighted 
in green (red) colors in tlie EDC panel. The EDC at the T momenta is in red color, c, In-plane 
spin-resolved EDC spectra, d, In-plane spin polarization obtained from c. 



Mn-Bi2Se3 film II is measured at the 13 beamline at Maxlab [Ml [35] in Lund, Sweden, 
with the spin-resolved energy distribution curve (EDC) mode (along binding energy, 
with fixed momenta). The in-plane and out-of-plane measurements are shown in Fig.js] 
and Fig. |4] respectively. Since manipulator of the 13 beamline is equipped with an 
additional motor to tilt the sample, the small tilt due to sample mounting is corrected 
before measurements. As a result P^, is zero in the measurements. Again, P^^ shows the 
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Supplementary Figure 4: Out-of-plane spin measurements of Mn-Bi2Se3 film II. In- 
plane spin measurements of Mn-Bi2Se3 film II a-b, Spin-integrated ARPES measured dis- 
persion and EDCs with incident photon energy 9 eV (used for spin-resolved measurements). The 
EDCs selected for spin-resolved measurements are highlighted in green (red) colors in the EDC 
panel. The EDC at the T momenta is in red color, c, Out-of-plane spin-resolved EDC spectra, d, 
Out-of-plane spin polarization obtained from c. 

"classical" helical spin texture. Moreover, shows the singly degenerate Dirac point 
and the TR breaking evidence, which has been presented in the maintext already. 

2. Surface magnetic interaction strength bz based on spin-resolved measure- 
ments on Mn-Bi2Se3 

Here we show the details of the method of fitting the magnetic interaction strength bz 
using the spin-resolved data sets on Mn-Bi2Se3, as presented in maintext Fig. 5c. 

We start from the simplest k-p hamiltonian for the topological surface states under 
TR breaking: 

fb vk \ 
^ +o{k'^) (1) 
vk_ -bz / 
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Supplementary Figure 5: Surface magnetic interaction strength based on spin- 
resolved measurements. The magnetic interaction strength can be fitted based on the spin- 
resolved measurements on films II and I (see maintext Fig. 3) of Mn(2.5%)-Bi2Se3. 

where bz is the magnetic interaction strength, /c+ and are defined as kj^ = ky + ik^ 
and k_ = ky — ik^ respectively. The higher k orders o(A;^) correspond to the parabolic 
bending of the linear dispersion at large momentum (/c^ term) and Fermi surface 
warping also at large momentum [k^ term), which are not significant for surface states 
at the energy near the Dirac point. This description is in particularly appropriate for 
Bi2Se3 since the surface Dirac cone is nearly isotropic and the warping effect is not 
strong as comparing to Bi2Te3. We solve the Schrodinger equation 



H|^ >= E\ij > (2) 
The eigen-values of the energy spectrum are 



E± = ±v/(t;A;)2 + (6,)2 (3) 

which represent the upper (E+) and the lower (E_) Dirac bands respectively. The 
eigen-states wavefunction for upper {\ip+ >) and lower {\ip_ >) Dirac bands are 



(4) 



1^- >= 



'^cos{9/2)\- 



> -e-^^/2^m(e/2)| 



> 



(5) 
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where |±| > are the eigen-states of electron spin in z direction, and 0) is parame- 
terized as 

{vkx,vky,bz) = a/ {vky + {bz)'^{cos{9)cos{(j)) , cos{9)sin{(j)), sin{6)) (6) 

The angles (6*, 0) exactly depict the spin texture, which are measured by spin-resolved 
ARPES. Where 9 is the out-of-plane polar angle of the spin vector, and is the 
in-plane azimuthal angle of the spin vector. From equation (S6), we obtain that 

9 = arctan(^) (7) 
vk 

This equation (S7) reflects the competition between TR-breaking effect and spin helical 
texture. We can use equation (S7) to fit the magnetic interaction strength bz as 
presented in Fig. 5c of maintext. 

As shown in Fig. |5| we have the fitting of bz using spin-resolved date sets from films 
II and I of Mn(2.5%)-Bi2Se3 respectively. The spin-resolved data on films I and II of 
Mn(2.5%)-Bi2Se3 are systematically presented in Fig. 3 of maintext. Fig. [2| Fig. |3| 
and Fig. |4] in the supplementary. For each measured spin vector, both in-plane spin 
polarization P// and out-of-plane spin polarization P^ are measured. The spin vector 
out-of-plane polar angle 9 is then experimentally obtained as tan9 = Thus for 
each experimentally measured spin vector, the data can be presented in the way of 
{9, k), where 9 is the experimentally measured out-of-plane polar angle 9, and k is the 
momentum location of that spin vector in the momentum space. Using these data 
sets, the magnetic interaction strength bz can be fitted by equation (S7). 

Fig. |5] shows the fitting based on spin-resolved date sets from films II and I of 
Mn(2.5%)-Bi2Se3 respectively. The fitted bz value is found to be 21 ± 5 meV for 
data sets fitting from film II, and 15 ± 6 meV for data sets fitting from film I. 

3. Remanent magnetization and temperature dependence of the TR breaking 
out-of-plane spin texture 

The spin-resolved data on Mn-Bi2Se3 shown in the main paper is done in remanent 
magnetization mode: the film is pre-magnetized in the UHV preparation chamber in 
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Supplementary Figure 6: Remanent magnetization and temperature dependence of 
the out-of-plane spin polarization of Mn-Bi2Se3 film surface, a, Out-of-plane spin polar- 
ization Pz with film magnetization in -z direction at T=45 K. Such condition is achieved by apply 
0.4T magnetic field in -z direction in the prep chamber and then transfer the film into the mea- 
surement chamber, b, Out-of-plane spin polarization with film magnetization in +z direction 
at T=45 K. Such condition is achieved by apply 0.4T magnetic field in +z direction in the prep 
chamber and then transfer the film into the measurement chamber, c, Out-of-plane spin polariza- 
tion at T=300 K. No magnetization is induced because temperature is above Curie temperature 
of the film surface (see below ferromagnetic order section). 

the ARPES system at the base temperature by approaching the film with a perma- 
nent magnet with magnetic field roughly at 0.4T, in order to induce the remanent 
magnetization. The the film is transferred from the preparation chamber to the mea- 
surement chamber to perform spin-resolved measurements, maintaining the pressure 
and temperature. 

Here we show the spin resolved measurements on Mn-Bi2Se3 film II as a function of 
remanent magnetization directions and temperatures. The helical in-plane component 
of the spin polarization does not show any dependence with magnetization or temper- 
ature. Thus we focus on the out-of-plane component which is a direct consequence 
of the magnetization of the film surface. To capture the essence, we only show the 
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out-of-plane spin polarization measurements at the TR breaking invariant T momenta 
k =0.0 A ^ . Fig. [g^ repeats the measurements shown in the maintext and in Fig. [ijl. 
When the remanent magnetization is reversed, as shown in Fig. [6]d, the out-of-plane 
spin polarization is also observed to be reversed. Identical measurements at room 
temperature 300 K (Fig. ^) show no observable out-of-plane spin polarization. 

4. First principles electronic structure calculation of Mn-Bi2Se3 and Mn im- 
purity band 



a Undoped BijSe^ b 4% Mn-doped Bi^Sej 




Supplementary Figure 7: First principles electronic structures calculation of undoped 
and Mn-doped Bi2Se3. a, undoped Bi2Se3. b 4% Mn-doped Bi2Se3. The band character of Mn 
Id and 2d orbital are plotted on the band structure calculation. 

Here we show our first principles electronic structure calculation of the Mn-doped 
Bi2Se3. Fig. [T] shows calculated band structure of both undoped and 4% Mn-doped 
Bi2Se3. Comparing to the undoped material, Mn doping is found to hole-dope the 
system, which moves the Fermi level (energy=0) below the Dirac point. The hole 
doping is due to the 2"*" ionic state for Mn atom (as compared to a S'^ ionic state for 
the Bi atoms in Bi2Se3). Such p-type doing behavior is also directly observed in our 
ARPES measurements. 

Apart from moving the chemical potential, an important feature in the calculation 
is the development of extra electronic states with strong Mn orbital character (Mn 
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impurity bands). The impurity bands are observed to be majorly located at energies 
lower than 0.3 eV below the Dirac poin in our calculation (see Fig. [7j3). At the same 
time, the valence band spectral weight is found to be shifted from ± 0.1 A 
to A;//~0.0 A ^ by Mn doping, as shown in Fig. Ic in the maintext, measured with 
the same photon energy. The newly developed strong intensity at the valence band of 
Mn-Bi2Se3 in ARPES measurements can have strong correlations to the Mn impurity 
bands in the calculation. 



5. Ferromagnetic ordering in Mn-Bi2Se3 film 
5.1. SQUID measurements of Mn-Bi2Se3 film 
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Supplementary Figure 8: Magnetic property of Mn-Bi2Se3 measured by SQUID, a, 

Mn-Bi2Se3 film sample layout. The indium layer on the bottom which is necessary for the sample 
mounting for MBE growth becomes superconducting below 4 K which contributes a very large 
diamagnetic signal to the SQUID measurements . b, Remanent magnetization measurements of 
Mn-Bi2Se3 film with zero field cooling (ZFC). The big drop of magnetization signal comes from 
the indium layer at the back side of the sample, c, Hysteresis measurements at T=1.5 K. The 
anomaly at H ~ again comes from the diamagnetic signal of the superconducting indium layer 
on the back side of the sample. 

As a starting point, the magnetic property of the Mn-Bi2Se3 film can be measured by 
superconducting quantum interference device (SQUID). As shown by the remanent 
magnetization measurements on the samples in Fig. [8]d, clear remanent signal can be 
observed for temperatures below 10 K. Clear hysteresis is also shown by SQUID below 
5 K (Fig. Isfc). The anomalous drop of the magnetization signal in both remanence 
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and hysteresis comes from the diamagnetic signal of the superconducting indium layer 
(below 4 K) on the back of the film sample, which is not relevant to the magnetic 
property of the Mn-Bi2Se3 film itself at all. 

In general, SQUID measurements are not so sensitive to small magnetization signals. 
In fact, small remanent signal is observed by SQUID above 10 K, but the signal is 
too small and no clear hysteresis is observed as a result. Moreover, SQUID includes 
magnetization signal of the entire sample (even the irrelevant indium superconducting 
diamagnetic signal on the back side of the sample). Thus a surface sensitive magne- 
tometer with a better sensitivity to small magnetization signal is needed to measure 
the magnetic properties on the top surface of the Mn-Bi2Se3 film, where the topologi- 
cal surface electrons are localized, as well as where ARPES and spin-resolved ARPES 
mainly probe (within ~ 20 A). 

5.2. XMCD measurements on the top surface of the Mn-Bi2Se3 film 

The magnetic properties of the surface of the Mn-Bi2Se3 film need to be treated care- 
fully for two reasons: (i) It has been theoretically predicted that the ferromagnetic or- 
der of the TI surface can be achieved via the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction mediated through the topological surface Dirac fermions jl2], even in ab- 
sence of the bulk magnetic ordering [15J. (ii) More importantly, careful Mn concen- 
tration characterization measurements have found that the surface has much higher 
Mn concentration than the entire film crystal |13] . The Mn concentration of the entire 
film crystal, which is the nominal concentration used in the maintext, is measured 
by the in-situ Bi/Mn fiux ratio during the MBE growth in conjunction with ex-situ 
Rutherford back scattering, both of which reveal the ratio of Mn atoms versus the 
Bi atoms over the entire crystal. On the other hand, the depth profile of Mn doping 
concentration can be measured by secondary ion mass spectrometry (SIMS), which 
shows that an excess of Mn accumulates on the surface during the sample growth. 
SIMS shows that for the top 50 A of the surface, the Mn concentration is about 8 
times higher than deep in the bulk. Finally, scanning tunneling microscopy (STM), 
which solely probes the top surface, also reveals a much higher Mn concentration on 
the surface as compared to the nominal Mn concentration, as well as the the fact that 
Mn atoms are inhomogeneously distributed on the surface, including the presence of 
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clusters of Mn atoms on the sample surface |13]. The high Mn concentration on the 
surface is very likely to be due to the diffusion of Mn atoms toward the film surface 
at high temperatures (> 200 K) during MBE growth or during decapping procedure 
used for ARPES and STM measurements. In the well-known dilute magnetic semi- 
conductors (Ga, Mn)As, the Curie temperature is found to be proportional to the Mn 
concentration. Mean- field Zener model Plj predicts TcOcx-p^/^, where x is the dop- 
ing concentration of Mn to Ga, and p is the hole carrier density. Therefore the high 
Mn concentration on Bi2Se3 film surface and its p-type doping to the Bi2Se3 system 
are both favorable to the ferromagnetic order on the surface at the relatively high 
temperatures(> 10 K). 

Here we use the total electron yield (TEY) mode of X-ray magnetic circular dichroism 
(XMCD, for review of the technique and its application to transition metal thin films 
and monolayers, see |l5l |36] ) to directly measure the magnetization of the surface. 

5.2.1. XMCD on Mn-BiaSeg film top surface at T=45 K 
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Supplementary Figure 9: XMCD measurements at Mn L23 absorption edge on Mn- 
Bi2Se3 film top surface at T=45 K. a, Remanent normalized XAS intensity of Mn L23 absorp- 
tion edge. Black (red) represents the XAS with the photon spin parallel (anti-parallel) with the 
remanent magnetization, b, The difference of the two XAS (parallel and anti-parallel) generates 
the XMCD signal, which is proportional to the magnetization of the Mn-Bi2Se3 film top surface, 
c, Hysteresis measurements are performed with field applied along out-of-plane direction. The 
measurements are done by measuring XMCD signal at the L3 edge energy E~638.2 eV at different 
H fields. 
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Fig. [9] shows the XMCD measurements on Mn-Bi2Se3 film top surface. The remanent 
measurements are performed at T=45 K and H=OT after applying magnetic field 
H=0.5T along -z direction. If the surface is ferromagnetically ordered at 45 K, then 
a finite remanent magnetization M should take place along -z direction after applying 
the magnetic field. Such remanent magnetization then will lead to a different response 
of the XAS spectrum using circular plus or circular minus polarized light, resulting in 
a nonzero remanent XMCD. Indeed, a clear difference between the different photon 
spin XAS spectra are observed, as shown in Fig. |9^. The L3 edge is observed to 
be enhanced with photon spin parallel to the magnetization, whereas the L2 edge is 
enhanced with photon spin anti-parallel to the remanent magnetization. The clear 
remanent XMCD (Fig. |9}d) suggests that the surface is ordered with easy axis along 
the out-of-plane direction at T=45 K. 

The hysteresis then can be measured by measuring XMCD signal (we choose the peak 
energy ofs L3 edge E~638.2 eV) at different applied H fields. As shown in Fig. a 
clear hysteresis is observed at T=45 K along out-of-plane direction, which proves the 
ferromagnetic order of the Mn-Bi2Se3 film top surface. 

5.2.2. XMCD temperature dependence 
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Supplementary Figure 10: Temperature dependence of remanent XMCD. XMCD mea- 
surements are performed at three different temperatures: T=45 K, 110 K and 300 K (room tem- 
perature) . 

Fig. [To] shows the temperature dependence of the XMCD signal. The XMCD signal is 
proportional to the remanent magnetization of the film surface, which can be obtained 
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by the sum rules [? ] (see below). Thus, the temperature dependence of the remanent 
XMCD signal reveals the temperature dependence of the magnetic ordering of the 



film surface. As shown in Fig. 10, clear XMCD is observed at T=45 K, whereas the 
XMCD is strongly suppressed at T=110 K. Finally no XMCD is observed from the 
measurement performed at room temperature. The detailed Curie temperature of the 
surface requires future finer measurements with more intermediate temperatures. 
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SI IV. Zn-BisSeg 



1. Detailed spin-resolved measurements on Zn-Bi2Se3 
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Supplementary Figure 11: In-plane spin measurements of Zn-Bi2Se3. a-b, Spin- 
integrated ARPES measured dispersion and EDCs with incident photon energy 9 eV (used for 
spin-resolved measurements). The EDCs selected for spin-resolved measurements are highlighted 
in green (red) colors in the EDC panel. The EDC at the f momenta is in red color, c, In-plane 
spin-resolved EDC spectra, d, In-plane spin polarization obtained from c. 

Here we show comprehensive spin-resolved measurements on Zn-Bi2Se3. These mea- 
surements on Zn-Bi2Se3 are performed at the 13 beamline at Maxlab [SllES] in Lund, 
Sweden, with incident photon energy 9 eV as a control group to the Mn-Bi2Se3. For 



in-plane spin polarization (see Fig. 11 ), the spin-helical configuration is again observed 



without surprise. On the other hand, the out-of-plane spin polarization (see Fig. 12) 
is zero at small momentum near the T momenta which is in sharp contrast to the Mn- 
Bi2Se3 case. At large momentum (k ~ —0.10 A) a very small but finite appears. 
However, the in Zn-Bi2Se3 respects TR symmetry since it reverses its sign upon 
going from —k to +k. As shown by Fig. [12)3, the Fermi surface of Zn-Bi2Se3 is warped 
into hexagon- shape. Therefore, the observed Pz should be a result of the Fermi surface 
warping on the Zn-Bi2Se3 Dirac cone [47]. 
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Supplementary Figure 12: Out-of-plane spin measurements of Zn-Bi2Se3. a, Spin- 
integrated ARPES measured dispersion witli incident photon energy 9 eV (used for spin-resolved 
measurements), b, Spin-integrated Fermi surface mapping of Zn-Bi2Se3 at the native chemical po- 
tential (Eb = 0.0 eV) c, Out-of-plane spin-resolved EDO spectra, d, Out-of-plane spin polarization 
obtained from c. 

SI V. Ultra-thin undoped Bi2Se3 

1. SR-APRES measurements, 3 QL vs 60 QL 

Here we show the spin polarization measurements on 3 QL undoped Bi2Se3 (Fig. 
4e-h in the maintext) again, and comparing the results with the spin polarization 
measurements on 60 QL undoped Bi2Se3. The thickness of 60 QL is well-above the 
quantum tunneling thickness limit. Thus no tunneling or coupling effect between the 
top and bottom surfaces is expected. This can also be seen from the ARPES measured 



gapless dispersion, as shown in Fig. 13:. The 3 QL and 60 QL SR measurements 



are performed with identical experimental conditions, except the film thickness. As 



shown in Fig. 13 i, no reduction of the spin polarization is observed on 60 QL film 



surface, even at small momentum such as ±0.03 A ^. Such measurements demonstrate 
the momentum resolution of the instruments, and show that the reduction of spin 
polarization on 3 QL film surface is not caused by blurring between the two branches 
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Supplementary Figure 13: SR-APRES measurements, 3 QL vs 60 QL, demonstrating 
the momentum resolution of the instrument, a-b, SR-ARPES measurements {hu = 60 eV) 
on 3 QL undoped Bi2Se3 shows reduction of spin polarization due at small momenta near the gap. 
Out-of-plane spin polarization does not show any significant spin polarization P^, especially at the 
vicinity of the f momenta, c-d, SR-ARPES measurements on 60 QL undoped Bi2Se3 measured 
with identical experimental settings {hv = 60 eV) does not show any reduction of spin polarization, 
e, in 3 QL ultra-thin film, the surfaces from top and bottom are quantum mechanically coupled to 
each other, resulting in a tunneling gap as observed in panel a. 



of the upper Dirac cone due to the momentum resolution of the instrument. 



